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Abstract 

Nonlinear properties of the order parameter modulation wave in such sys- 
tems as thiourea are described in the framework of the phenomenological 
model with no Lifshitz invariant. It is also shown that for some values of the 
thermodynamic potential parameters the lock-in transition is continuous in 
the type II systems. Close to this transition the soliton lattice appears. 
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1 Introduction 

When describing the incommensurate (IC) phases in the ferroelectrics of so called 
type II im the phenomenological approach proposed by Y Ishibashi and H Shiba 
is widely used 0. 

In according to this model, the system thermodynamic potential can be written 
in the form |^: 

L 



where '^[x) is a one- dimensional order parameter (e.g., a component Py[x) of the 
spontaneous polarization P =Py); = dip/dx, L is a crystal length in the direc- 

tion of spatial modulation of the order parameter. 



remain the notation '7' to indicate the contribution of invariant (v^')^ which favors 
the appearance of IC state). For all known ferroelectrics of type II (sodium nitrite 
NaN02 thiourea SC{N 112)2 05 Sn2P2SeQ and betaine calcium chloride 
dihydrate (BCCD) 0) the parameter g is negative: g < 0. For sodium nitrite 
and Sn2P2SeQ the direct (virtual) disordered-to-commensurate phase transition is 
regarded to be of the first order, thereby p < and = 1 for NaN02 P and 
Sn2P2SeQ 0. In the case of thiourea and BCCD the parameter p > and h = 0. 
It is usually assumed that the only parameter q is dependent on temperature T: 
q = qo ■ (T — To), where go and Tq are some constants. 

The model (|1]) fairly well describes a lot of properties of the IC phase in the 
ferroelectrics of type II 0. One of such properties is a predominantly sinusoidal 
character of the order parameter modulation wave. Experimental studies indicate 
that for the type II ferroelectrics higher order satellites are of low intensity even in 
the close vicinity of the lock- in transition (see for review 0). This circumstance is 
a reason why the one-harmonic approximation is often used when considering the 
IC order parameter configuration in, e.g., sodium nitrite or Sn2P2SeQ [0, H, H, 

ip{x) = a ■ sm{bx). (2) 

On the other hand, the recent experiments |T^, reveal that for thiourea the 
dependence of order parameter on position x contains relatively large contribution 
of the higher harmonics and can not be regarded as pure sinusoidal. 

The solitonic properties of the modulation functions in thiourea are explained in 



1^ . In order to describe nonlinear features of the IC order parameter configuration 



the authors of make use another than (|1]) phenomenological approach |T3 

The approach [1^ is similar to the model developed for the systems with two- 



component order parameter for which the Lifshitz invariant can be introduced |n 
Examples of such systems are compounds of the A2BX4 family, e.g., Rb2ZnCU 
[y]. Although in the case of type II ferroelectrics the order parameter is usually 
considered as one-dimensional, the role of second component is performed by some 
other normal coordinate ^(x) (e.g., xy-component of the elastic strain tensor) |T3[| . 
The function ^ (x) transforms like the first derivative of the order parameter: ^ (x) ~ 
(p'{x). These transformation properties of ^(x) allow to construct the thermodynamic 
potential with a term analogous to Lifshitz invariant [|13| . 



The approach [|l^] is expected to be more appropriate than the model (0) when 



interpreting nonsinusoidal configurations of the order parameter |T^. To some ex- 
tent, these expectations are grounded on the analogy between the approach [|13 



and the theory developed in \T^\. Really, the latter constitutes a powerful tool for 
description of the soliton structures in the compounds of A2BX4 family |T^. 

In the present paper we show that the nonlinear contribution to the modulation 
wave observed in the experiments for thiourea can be explained in the framework of 
the model (|lD as well. 

Our consideration is based on the nonlinear approximation for the IC order 
parameter configuration proposed in [Q. Using it we obtain that if the system 
material parameters g and p are so that g ■ p~^ ~ —6.0 {h = 0, p > 0) then the 
ratio Os/ai for the amplitudes of third (03) and fundamental (ai) harmonics of the 



Moreover, it is found that for some values of the parameters g and p (for systems 
with = 0, p > these values satisfy the relationship ■ gp~^ ~ —16) the transition 
from the IC phase into the commensurate state is continuous. The order parameter 
configuration is domain-like in the proximity of such transition. 

The structure of present paper is following. In section 2 we formulate main 
features of the order parameter approximation proposed in [Q. In section 3 we 
consider the nonlinear configurations of order parameter in the case of thiourea. 
The estimations of nonlinear properties for other known ferroelectrics of type II are 
given as well. The principal possibility of existence of the strong soliton regime in 
the type II systems is investigated in section 4. The comparison of the obtained 
results with the properties of type I systems is presented in section 5. 

2 Sn-approximation for the IC order parameter 
configuration 

In according with the analysis made in |Q, the equilibrium configuration of order 
parameter in the IC phase for the type II systems (|^) can be approximated as 

(p{x) = a ■ sn{bx, k) (3) 

where sn{x,k) is the Jacobi elliptic sinus |jl6[. In (^) the amplitude a, the wave 
number b and the elliptic modulus k {0 < k < 1) are defined by the minimization 
of the thermodynamic potential (|I]) in respect to a, 6, 0]. 

In contrast to the approach (|^), the approximation (^ allows to consider not 
only the linear regime of the IC phase but the nonlinear configurations of order 
parameter as well. 

Really, if the elliptic modulus k is small {k ^ 0) the function (|^) is closed to the 
dependence (H): sn(x. A; 0) ~ sin(a;) [|l^. But when k ^ 1 the spatial behavior 
of elliptic sinus becomes domain-like: the wide regions with almost constant values 
(p{x) ~ ±ipo are separated by the narrow regions where the function (^ changes 
abruptly. 

In the model (|1|), (^ the elliptic modulus k is equal to zero at the point qj = 
of the disordered-to-incommensurate phase transition 0]. With decreasing tem- 
perature the elliptic modulus k grows and takes its largest value kc at the lock-in 
transition point Qc- The preliminary investigations show that k^ can be close to 
unity for some values of the material parameters. For example, ii g = —10, p = 1, 
h = then k^ = 0.965 g]. 

Another important property of the order parameter approximation (^ is an 
additional (in comparison with the approach (0)) mechanism causing the change of 
modulation period. 

The continuous dependence of modulation period P on temperature is one of the 
most characteristical features of the IC phases 0. In the framework of approach 
(0) the period of IC structure is equal to P = 2ti /h. The wave number h depends 
on temperature only when the {^p^p')"^ - invariant is present in the expansion of the 
thermodynamic potential (|I|) |^ 



And if the material parameter g is negative then the period P increases with de- 
creasing temperature as it is observed in experiments ||^. 

The period of IC order parameter configuration (j^) is defined as P = 4K{k)/b 
{K{k) is the complete elliptic integral of the first kind |jl6[). It depends not only on 



the value of h (as it takes place in the one-harmonic approach @)) but also on the 
elliptic modulus k which characterizes the degree of nonlinearity of the modulation 
wave. As consequence, the approximation imposes less strong requirements on 
the material parameters, in particular, on g (e.g., the period P grows even when 
(7 = 0, Moreover, if the elliptic modulus k is close to unity k 1 (i.e. 

K{k) oo [|T6|) then the nonlinear mechanism of the increase of modulation period 
becomes dominant and P can be very large: P ^ oo. 

The numerical investigation of variational equation for the functional (|T]) 
shows that the approximation (|^) correctly reproduces nonlinear properties of the 
equilibrium configuration of order parameter in the IC phase. 



3 Nonlinear configurations of the order parame- 
ter in thiourea 

Now let us to apply the model (|Tp, (§P to thiourea. 

As it has been mentioned above, in the case of thiourea the material parameter 
h equals to zero. For the sake of simplicity we also assume that p = 1. 

Numerical minimization of the thermodynamic potential (|I]), (0) in respect to a, 
b, k shows that if the material parameter g is equal to g = —6.0 then the elliptic 
modulus takes the value kc ~ 0.923 at the point Qc = —0.290 of the lock-in transition. 
Using the formulae for the Fourier expansion of elliptic sinus or the procedures 
of fast Fourier transformation (FFT), one can find that for this kc the ratio of 
amplitudes of the third (03) and the fundamental (ai) harmonics of the modulation 
wave is equal to a^/ai ~ 0.104 and also a^/ai ~ 0.012 (05 is an amplitude of 
the fifth harmonic). The spatial dependence of the order parameter (p{x) at the 
temperature Qc is shown in figure 1 (full line). The obtained results are in a good 
agreement with the experimental data |iy, |lT|]. 

Besides the order parameter f{x), some other modulation functions are also 
discussed in the case of thiourea |]T2[. If such function ^(x) is coupled with the 
order parameter ip{x) = asn{bx,k) = asin[6'(x)] by the relation ^(x) ~ cos[6'(a;)] = 
cn{bx, k) then the dependence of C,{x) on position x has the form depicted in figure 
1 by broken line (cf. with figure 3 in [0). 

Of course, the complete description of ^{x) and other modulation functions re- 
quires a modification of the thermodynamic potential (|I]) by including additional 
invariants which correspond to interactions of C,{x), ... with each other and with the 
order parameter. Such an analysis is, however, beyond the scope of the present 
paper. 

The nonlinear properties of the IC modulation wave in thiourea essentially differ 
from ones for sodium nitrite. Taking for estimations the material parameters given in 
P (in our notations they correspond to g = —9.51, p = —0.651, h = 1) we find that 



as/ai ~ 0.0007. Note that these resuhs are very close to ones obtained analytically 
in 



19 



As for other ferroelectrics of type II, estimations reveal that for Sn2P2SeQ {g = 
—1.37, p = —0.19, h = 1 0) the characteristics of modulation wave at the lock- 
in transition point are the following: kc ~ 0.706, a^/ai ^ 0.041, a^/ai ~ 0.0018; 
and for BCCD {g = -8.0, p = 2.0, h = ^) - ^ 0.887, a^/ai ^ 0.086 and 
as/ai ^ 0.0082. 



4 Soliton regime in the type II systems 



Results given in the preceding section describe the nonlinear properties of the IC 
modulation in four known compounds belonging to the type II ferroelectrics. How- 
ever, these estimations do not answer on more general question: to what extent the 
soliton regime can develop in the model (|I|), in principle. In order to clarify the 
situation we have investigated the behavior of systems (|) in the limit k ^ 1 in 
more detail. 

When k ^ 1 the thermodynamic potential (|l|), (|^) acquires the form [Q: 



\p (1 + k'^) - A-i (|p + yh^) - k'^K-^ [Ip + ^ 
|/..(l + |fc--iA-i-ifc-A-) 



+ 



(5) 



+ 



where k''^ = I - k\ A= ln(4/A;'); A^^ -> if A; ^ 1. 

The equilibrium wave number b can be found from the equation d^/db 
what yields (cf. with (|)): 



-7 ■ (1 + - 

8 ' ^ 2 



k'^) + -g-a^ 



Taking into account 



1 + k") 

the thermodynamic potential (Bl) 



(6) 



can be rewritten as 



1 



IU2 



q - A-i 
ip-A-i (Ip+j-^ig 



+ 



(7) 



(1 + k'^ 

{l + Ik-y. [Ih - A-^ {§h + ^,g 



The function A ^{k) changes much slower then k'"^ approaches to zero (e.g., for 
k'^ = 0.1 the value of A-\k) is 0.394, for k'^ = A x 10"^ - A-\k) ^ 0.1). Hence, it 
is reasonable to omit the terms proportional to k'^ in (|^). 

The further analysis depends on the value of material parameter h. For the 
systems with h = 0, p > the results are as follows. 

The equation d^/da = defines the equilibrium amplitude a. Substituting its 
solution into the expression (|^) and comparing the result with the thermodynamic 
potential for the commensurate state $c = — 9^/2p, we find the effective temperature 
Qc of the lock-in transition: 



In (^) the approximate (right hand side) formula is derived for the case when the 
term proportional to is neglected in (^. This expression clearly demonstrates 
main features of the dependence of qc on gp~^. 

Now we consider the conditions under which the elliptic modulus can be equal 
to unity at the lock-in transition point Qc- 

The thermodynamic potential (0) depends on the elliptic modulus k only through 
the function A~^(fc) (remember we omit the terms proportional to fc'^ in (j^)). 
Thereby it is convenient to formulate the variational task for k in terms of A~^. In 
these terms the condition /cc = 1 means hr^{k^ = 0, and the equation d^/dk = 
is equivalent to d^/d{A~^) = 0. Solving the latter at the temperature qc (||), we 
find that the function A~^{kc) is equal to zero if the material parameters satisfy the 
relation: 

^.gp-' = -lQ. (9) 

Therefore, for the systems with material parameters g and p related with each 
other in according to (|^) {h = 0), the elliptic modulus of the order parameter 
modulation wave (^ is equal to unity at the point of lock-in transition. It means 
that for such systems the transition from the IC phase into the commensurate state 
is continuous. Close to this transition the order parameter spatial configuration is 
domain-like, and the soliton density ns = 7^/2K{k) (see [1^ and references therein) 



approaches to zero with decreasing temperature q to its lock-in value qc- 

Numerical calculations confirm the results of analytical investigation (see table 

!)• 

For the case gp ^ = — 16 the elliptic modulus kc has approached to unity, but 
we have stopped calculations at the value kc = 0.9999. The properties of elliptic 



functions change abruptly at the point k = 1 , and there are some difficulties to 
reproduce the point k = 1 numerically. 

As it follows from table 1 the dependence of kc on gp^^ has a maximum when 
gp"^ = —16. For little {gp~^ ~ 1) and large {gp^^ ~ 100) values of the parameter 
combination gp^^ the contribution of higher harmonics in the modulation wave is 
relatively small. 

The analysis for the case of systems with h = 1, p < can be made in the similar 
manner. Here we point out only the following. 

When h = 1, p < 0, the direct disordered-to-commensurate phase transition is 
of the first order. As consequence, for large enough values of the parameter \p\ the 
range of IC phase stability can be relatively small like it takes place for NaNO^ 
[qc ~ +0.07) 1^. Nevertheless, for any p there exists some g for which kc = 1. For 
example, for p = —0.65 (the case of sodium nitrite 0, ^) the lock-in value kc of the 
elliptic modulus of the modulation wave (^ equals to 1 if the material parameter g 
is g ^ —4.5. 

However, due to the presence of term proportional to in (|I]) , i.e. due to 
h = 1, the material parameters g and p are not so correlated as in the case of h = 
(in fact, when h = the ratio g/p is relevant rather than the parameters g and p 
themselves). As consequence, for the systems with h = 1, p < the dependence 
g{p) providing kc = 1 is more complex than (^. 



5 Discussion 



In the present paper we have shown that in the framework of the phenomenological 
model with no Lifshitz invariant |^ different nonhnear configurations of the IC 
order parameter can be described: almost sinusoidal one as in the case of sodium 
nitrite (03/01 < 0.03); one with more large contribution of higher harmonics as 
in thiourea (03/01 ~ 0.1); strong soliton regimes when the lock-in transition is 
continuous (such compound is unknown at the moment). Nonlinear properties of 
the concrete system are defined by values of the material parameters g and p (see, 
e.g., table 1). 

The specific role of {(pf')"^ - invariant should be emphasized when discussing 
the nonlinear features of the IC order parameter configurations. If this term is 
not included in the thermodynamic potential (|l]) the soliton structure does not 
develop (see also 0, |17|). From this point of view, (v?V^')^ " term is analogous to the 
Umklapp invariant of relatively low order (the anisotropy invariant) which favors 
the appearance of domain-like structures in the type I ferroelectrics [0, The 



difference between (</?</?') ^ - term in ([^) and the anisotropy invariant is that the 
former is a part of the gradient gain of the thermodynamic potential and can not 
influence on characteristics of the commensurate phase. On the contrary, in the type 
I ferroelectrics the behavior of IC and commensurate phases is correlated due to the 
anisotropy invariant (it belongs to the local interactions). In the case of type II 
systems analogous correlation appears only for some specific values of the material 
parameters g and p, i.e. only for some correlated actions of the invariants {ff')'^ 
and which define the low-temperature behavior of IC and commensurate phases. 

The existence of IC state in the type I systems is caused by symmetric prop- 
erties (the Lifshitz condition is not fulfilled) 0]. For the systems of type II such 
global reasons are absent and the spatial modulation of the order parameter is a 
consequence of specific features of the interatomic interactions . 

As result, in the type II ferroelectrics the appearance of soliton regime has no 
systematic character, in contrast to the situation which takes place for, e.g., com- 
pounds of the A2BX4 family. 
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Figure Captions. 



Figure 1. The order parameter (p{x) (full line) and the modulation function 

^{x) (broken line) as functions of position x at the temperature Qc — —0.29 (the 
lock-in transition point). The material parameters are the following: g = —6.0, 
p = 1.0, h = 0. The amplitudes of order parameter (p{x) and function ^{x) are 
arbitrary. 



Table 1. Characteristics of the order parameter modulation wave in the type 
II systems (the case /i = 0, p > 0) at the point Qc of lock-in transition for different 

values of the material parameters g and p: kc is a lock-in value of the elliptic modulus, 
tts/ai is a ratio of the third and first harmonics, ns is a soliton density. 



Table Caption. 



g/p 



-1.0 
-10. 
-16. 
-20. 
-100, 



0.690 0.768 0.052 0.81 

0.212 0.965 0.136 0.57 

0.156 0.9999 0.262 0.27 

0.134 0.970 0.142 0.56 

0.048 0.766 0.052 0.81 
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